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ABSTRACT 

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel 
ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration 
measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, 
/3, the Balmer decrement, and the far infrared (IR) to 150 nm far ultraviolet (FUV) luminosity 
ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order 
to identify possible systematics in these various measures. We reiterate the finding of other 
authors that there is a large scatter between the Balmer decrement and the parameter, and 
that f3 may be poorly constrained when derived from only two broad passbands in the UV. 
We also emphasise that FUV derived SFRs, corrected for dust obscuration using /3, will be 
overestimated unless a modified relation between /3 and the attenuation factor is used. Even 
in the optimum case, the resulting SFRs have a significant scatter, well over an order of mag- 
nitude. While there is a stronger correlation between the IR to FUV luminosity ratio and /? 
parameter than with the Balmer decrement, neither of these correlations are particularly tight, 
and dust corrections based on /3 for high redshift galaxy SFRs must be treated with caution. 
We conclude with a description of the extent to which the different obscuration measures are 
consistent with each other as well as the effects of including other galactic properties on these 
correlations. 
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1 INTRODUCTION 

Dust obscuration in galaxies is a well-recognised and long stud- 
ied phenomenon. From the perspective of large galaxy surveys, the 
primary concern is often how to make suitable corrections for the 
dust obscuration within the galaxies of interest, in order to estab- 
lish their intrinsic properties such as SFR and stellar mass. This 
has traditionally been accomplished through the use of the Balmer 
decrement (Osterbrock 1989), a dust sensitive emission line ratio 
that is sU'aightforward to measure in optical spectra. With the ad- 
vent of ultraviolet (UV) and infrared (IR) satellite telescopes, new 
approaches to measuring or constraining global galaxy dust prop- 
erties have been established, with common metrics being the UV 
spectral slope, /3 (Me urer et al.,|1999| >, and the far infrared (FIR) to 
ultraviolet luminosity ratio, Lfir/Luv (Bell|2003) . 

The UV spectral slope f3 has been proposed as a suitable tool 
for deriving obscuration corrections, in particular for galaxies at 
high redshift where the Balmer decrement is not easily measurable 
(e.g., |Bouwens et al.,|20 09l, having been shifted to infrared wave- 
lengths. The effectiveness of /3 as an obscuration metric depends 
on how well the UV slope can be measured ( C alzetti et al.,|1994j 
|Kong et al.,|2004[ >. This is influenced by factors such as the instru- 
mentation used and associated sampling of the observed UV spec- 
trum, together with the source redshift impacting on the rest-frame 
UV wavelengths being probed. From recent comparisons of f3 with 
other obscuration measures, using samples selected at optical ( |Wi-| 
jesing he et al.,|201(ty or FIR ( |Buat et al.,|2010| > wavelengths, it is 
apparent that there are significant limitations in the use of the (3 
parameter for making obscuration corrections. 

The P parameter was formalised by Me urer et al.,| J 1999} 
by exploiting the relationship between the ratio of the far-infrared 
(FIR) and the ultraviolet (UV) fluxes and the UV spectral slope for 
a sample of 57 starburst galaxies. |Meurer et al.,|(T999| > argue that 
since the FIR flux in starburst galaxies is produced from the UV ra- 
diation that is absorbed and re-radiated by dust, the FIR to UV flux 
ratio can be used as a measurement of dust absorption. Subsequent 
work has emphasised that this relationship, originally calibrated for 
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starburst galaxies, does not work so well for the general population 
of star forming galaxies (e.g., Bell 2002, Kong et al., 2004, Buat 
|et al.,|20 05 1, with typical galaxies deviating from the relationship 
between j3 and Lfir/Luv established by |Meurer et al." 1 1999| >. 

This is to be expected, as in quiescent spiral (non-starburst) 
galaxies both the assumption that the dust emission is solely pow- 
ered by UV photons and the foreground screen assumption breaks 
down. Thus as much as 30% of the dust emission in spiral galaxies 
is predicted to come from the old stellar populations i Popescu et] 
|al.,|201l] >. Furthermore, as shown in Popesc uet al.,|(201 \\ , at least 
10% of the dust heating could come from the old stellar popula- 
tion in the bulge, for a typical spiral with a bulge-to-disk ratio of 
around 0.3. Radiative transfer modelling of the distribution of stars 
and dust has also shown that the finite exponential disk distribu- 
tion of stars and dust, as opposed to the foreground screen distri- 
bution, could lead to significant differences in the determination of 
the attenuation of stellar light (Tuffs et al., 2004, Pope scu & Tuffs,| 
|2004[ l. Not least, spiral galaxies have both a large scale distribution 
of stars and dust as well as a clumpy component associated with 
the star forming complexes in the spiral arms ( Sauvage et al.,|2003| 
[Popescu & Tuffs, 2010). This renders a proper calibration of the 
P relation to the radiative transfer models challenging. Neverthless 
more empirical approaches have been taken, in the hope to find fur- 
ther parameters for calibrating this relation. 

|Kong et al.,| f2004), therefore, point out that the relationship 
between the IR-UV ratio and the UV slope is dependent on the 
star formation history as well as the dust content, implying that a 
straightforward application of this relation might not be accurate. 
They present a modified relation which accounts for the star for- 
mation history as well as the dust content. Wijesinghe et al., (2010) 
and |Buat et appOlO) both find that SFRs derived from the FUV 
corrected for obscuration using /? and the attenuation relation of 
Meurer et aQ(l999) are significantly overestimated. Wijesi nghe et| 



al„ H2010F propose a revised attenuation relation for /3 which elim- 



inates this overestimation. Significant scatter, more than an order 
of magnitude, however, remains between the SFRs inferred using 
dust corrections based on j3 compared to those based on the Balmer 
decrement. 

While the use of /3 is attractive as a measure of dust absorp- 
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tion at high redshift, in practice there are several considerations 
that need to be taken into account. The main assumption made 
by |Meurer eTaLl jT999> was that high redshift galaxies will have 
similar multiwavelength properties to local starbursts. It is unclear 
how the intrinsic spectral slope will be affected with increasing red- 
shift, and indeed recent evidence for very blue UV spectral slopes 
at the highest redshifts JBouwens et al.,|2009||2010||Bunker et al.,| 
|2010[ > is now being interpreted as evidence potentially for either ex- 
tremely low metallicity or variations in stellar initial mass function. 
This opens the question of how much intrinsic variation there is in 
the underlying UV spectral slope amongst galaxies at progressively 
lower redshift, in particular if recent suggestions for variation of the 
stellar initial mass function (Dave 2008, Wilkins et al., 2008a b 
Meurer et al., 2009, Gunawardhana et al., 2011) are borne out. If 
there is indeed broad variation in the UV spectral slope, then there 
needs to be a redshift, or potentially a SFR or other dependence 
included in the use of /? when making obscuration corrections. 

To address these issues, we explore the relationships between 
the FIR luminosity, the UV luminosity, the B aimer decrement and 
ft, all as a function of galaxy mass, SFR, and redshift. In §[2] wc 
describe the data used in this analysis. §[5] presents an analysis of 
the ft and Balmer decrement relationships, and §[4] explores dif- 
ferent approaches to parameterising the UV attenuation. §|5] in- 
vestigates the relationships between the FIR and FUV luminosi- 
ties. These measures are discussed and analysed in §[7] and in §[5] 
we summarise our results and conclusions. Throughout, all magni- 
tudes are given in the AB system, and we assume a cosmology with 
H =70kms- 1 Mpc~ 1 , Q M = 0.3 and Q A = 0.7. 



2 DATA 

The Galaxy And Mass Assembly (GAMAj^Jsurvey is a multiwave- 
length imaging and spectroscopic survey covering ~ 144 square 
degrees of sky in three 12° x 4° regions <|Driver et al |[ 2009 



[20t0)|Robotham et al.,|2010||Baldry et al.,|2010[|Hill et al.,|2010| > 
GAMA provides the redshifts, emission line measurements and 
UV/optical/near IR (NIR) photometry used in this analysis, with 
spectroscopy from 2dF/AAOmega on the Anglo-Australian Tele- 
scope and imaging from GALEX, SDSS and UKIDSS. The FUV 
band in the GALEX filter has an effective wavelength of 1528 A 
and the NUV band has an effective wavelength of 2271 A. K- 
corrections were applied to the observed GALEX UV magnitudes 
using KCORRECT.V4.1.4 dBlanton et al.,||2003| to infer the rest- 
frame magnitude at 1528 A. The galaxies in this analysis have a 
GALEX FUV magnitude range of 17.2 < rriFuv < 26.8, and a 
redshift limit of z < 0.35 due to the requirement for Ha to be in 
the observable spectral range. Stellar masses are calculated for the 
GAMA galaxies using spectral energy distribution (SED) models 
based on the ugrizphotometry bands (Taylor et al.,|20f0) |. 

The Hersche^jATLAS 2 (H- ATLAS f\ survey will ultimately 
observe 550 square degrees of sky in 5 bands (100, 160, 250, 350, 
500 ^m) detecting ~ 200, 000 galaxies spanning < z < 4 (Eales 
et al. 2010). The H-ATLAS Science Demonstration Phase (SDP) 
field covers an area of ~ 4° x 4° centered at a = 09 h 05 m , 8 = 
00° 30' in the GAMA 9hr field. It was mapped in parallel mode 



using the PACS (Pogl itsch et al.,|2010) and SPIRE (Griffin et al. 
2010) instruments on-board the Herschel Space Observatory jPil-| 
|bratt et aT~j 2010l. The 5-band Herschel maps and catalogues for 
this SDP field are described in |Ibar et aT7~J ( |2010| >, |Pascale et al.,| 
\2Q\Q\ and |Rigby et al.,| ( |20l6"| >. Sources in the > 5a 250 fim cat- 
alogue were matched to the SDSS DR7 (Abazajian et al. 2010) 
and thence to GAMA using a likelihood ratio method pSmifh et ak~J 
l |2010) resulting in 2423 250 /im sources with optical identifica- 
tions from SDSS to r < 22.4 at > 80 percent reliability. Of these 
sources, 1050 were found to be in common with the GAMA survey 
in the ~ 12.7 square degrees of overlap. From the reliability values 
calculated for these objects, we expect to have a 1.9 percent false 
ID rate in this sample. 

The IR and FUV luminosities that appear independently have 
also followed the same treatment. The energy balance method of da 
Cunha, Chariot & Elbaz (2008) is used to derive the IR (dust) lumi- 
nosity from 3-1000/xm by fitting to the the UV-sub-mm spectral en- 
ergy distributions of the galaxies. A large, stochastic library of stel- 
lar and dust emission models was used that includes a wide range 
of star formation histories, metallicities, dust attenuation, dust tem- 
peratures and different contributions by various dust emission com- 
ponents to the total infrared emission (polycyclic aromatic hydro- 
carbons (PAHs), hot mid-infrared continuum, and dust in thermal 
equilibrium; full details are given in in |da Cunha et~aLl [20081 ). 
Throughout, we use the notation Lm to represent the IR luminos- 
ity described above. Figure [T] shows the UV and 250 pm fluxes, 
colour-coded by the sum of the FUV and IR luminosities, along 
with lines of constant flux ratio, to illustrate the range of values 
that this sample probes. The sum of the FUV and IR luminosities 
are derived from frequency independent luminosities of these two 
quantities (units of solar luminosities). 

The data used in the SED fitting includes UV photometry 
from GALEX (Seibert et al. in prep.), optical ugriz from SDSS 
DR7 (Abazajian et al. 2009), infrared YJHK from UKIDSS LAS 
(Lawrence et al. 2007) all of which is re-measured from the im- 
ages in matched apertures defined in the r-band by GAMA (Driver 
et al. 2010; Hill et al. 2010). The optical/NIR fluxes have been 
PSF (point spread function) matched l |Hill et aT. , 2010), and the 
UV fluxes were derived using the the technique of Robotham et| 
[al^pOTT] ). See Smith et al. (in prep.) for a full description of the 
implementation of this method to H-ATLAS data. 

Of the 1050 GAMA sources detected within the H-ATLAS 
SDP region, 899 have Ha emission, and 875 have FUV measure- 
ments from the GALEX-MIS survey. Of these galaxies 221 have 
measurements from both the PACS (a > 5) and SPIRE instru- 
ments on the Herschel telescope while the rest have only measure- 
ments from the SPIRE instrument. Within the GAMA survey there 
are 96231 galaxies with Ha measurements and 104681 with FUV 
measurements. Where we show measurements that do not rely on 
having IR data, we use the full GAMA sample for which Ha emis- 
sion and UV photometry have been measured, after excluding AGN 
based on spectral emission line ratio diagnostics (Baldwin et al., 
|2004| |Kewley et al. 2001 ). Where the IR luminosities are required, 
we are limited to showing only those within the H-ATLAS SDP 
region. 



1 http://www.gama-survey.org/ 

t Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA 
§ http://www.h-atlas.org/ 



3 BALMER DECREMENT AND UV SPECTRAL SLOPE 

The Balmer decrement is the ratio of the stellar absorption cor- 
rected Ha and H/3 fluxes. For the GAMA emission line measure- 
ments, the stellar absorption corrections are applied as: 
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Figure 1. 250/xm flux compared against the FUV flux, illustrating the flux 
limits of the GALEX and H-ATLAS measurements in the current analysis. 
Both quantities are in units of Wm — 2 , and the lines are constant values 
of log(F25o/i ? Fuv)- The colour indicates the total of the IR and FUV 
luminosity. 
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where f a and f s are the observed and stellar absorption corrected 
fluxes respectively, EW is the equivalent width of the emission line 
and EW C is the correction for stellar absorption, taken to be 0.7 A 
( |Gunawardhana et al.,||201 T| l. [Gunawardha na et al.,| (2011) show 
that varying EW C for each galaxy will not have a significant effect 
on the observed trends. They show that a 5% difference in Ha EWs 
are observed, and only for the extremely low (log[Ha EW] < 0.9) 
Ha EWs. For this same reason, varying EW C values for different 
emission lines will also have a minimal effect. We, therefore opt to 
use a standard EW C value of 0.7A. The Balmer decrement is then 



BD = f Ha /f, 



Hf), 



(2) 



where BD is the Balmer decrement and fua and fsg are me ste l~ 
lar absorption corrected fluxes of the Ha and H/3 emission lines. 
Aperture corrections were performed on emission line data follow- 
ing the prescription of Hopkins et al. (2003; see also Gunawardhana 
etal. 2011). 

All Balmer decrements below the Case-B value of 2.86 ( |Os 
terbroc k|l989) were set equal to 2.86 as suggested by |Kewley et 
al-lpOOoj . The application of the Balmer decrement in correcting 
the Ha luminosities for dust obscuration is described in Wijesinghe 
|et al.,| < [2Tjl"0| >, and we use the same method here. In particular we 
adopt the obscuration curve of |Fischera & D opita ( 2005 ), shown by 
| Wijesinghe et al.,| ( |2010| > to be the most effective at producing self- 
consistent SFR estimates simultaneously from FUV, NUV, [Oil] 
and Ha luminosities. 

The j3 parameter is calculated using the observed fluxes ob- 
tained through the GALEX FUV and NUV filters. The UV spectral 
slope is determined from a power-law fit to the UV continuum of 
the form: 



fx oc X p , 



(3) 



where fx is the flux density per wavelength interval and A is the 
central rest wavelength ( jMeurer et al.,fT9"99"] >. Following |Kong et] 
|ap (2004) we use: 



log f fuv - log f nuv 
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Figure 2. (a) Comparison of Balmer decrement and the f3 parameter. There 
is a very weak trend, although for any given value of /3 the Balmer decre- 
ment can vary by factors of typically 2-3, but ranging up to ten. The blue 
line shows the mean values for the Balmer decrement in bins of /3 while 
the green line shows the line of best fit. The error bars show the median er- 
rors, (b) These histograms show how the Balmer decrements are distributed 
across three /3 ranges. As we move to higher /3 values it is clear that the typ- 
ical Balmer decrement also increases. The 3 bins contain 2510, 9788 and 
1115 galaxies respectively in the order of lowest to highest /3 ranges. Two 
bins (-1.5 < f} < -0.5 and 0.5 < fl < 1.5) are omitted for clarity, 
although they follow the same trend. 



where A f u v = 1528 A and Xnuv = 2271 A are the effective 
wavelengths of the far and near UV filters of the GALEX satellite 
and }fuv and /nuv are the mean flux densities per unit wave- 
length through these filters. 

Figure [2] shows a very weak trend between the Balmer decre- 
ment and /3, but this is dominated by the remarkably large scatter 
in the distribution. The least-squares fit shown, with an unremark- 
able correlation coefficient of r = 0.18, is BD = 0.64,3 + 7.64. 
While we do not claim any formal correlation between these two 
measures, we use this result below in an exploration of different 
parameterisations of the FUV attenuation, Afuv- In order to min- 
imise any bias in the least squares fit any Balmer decrements with 
values ^2.86 were excluded from the fit. This removes ~ 21% of 
the sample for the purposes of this fit. 

While there is a large scatter, it can be seen from Figure [2] 
that at the highest f3 values there is a concentration of systems 
with higher Balmer decrements. In other words, systems with low 
Balmer decrements tend to have flatter UV spectral slopes. It is 
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also interesting to note that the lowest values of (3 (the bluest UV 
spectral slopes) are not associated with the lowest Balmer decre- 
ments. Galaxies with the highest measured Balmer decrements tend 
to have a broad range of UV spectral slopes. Systems with flat UV 
slopes (/? ~ 0), however, have the broadest range of Balmer decre- 
ment. 

The Balmer decrement and the /3 parameter both measure the 
strength of obscuration in different ways, so naively a correlation 
between the two parameters might be expected. The dust geometry 
in galaxies is not a simple foreground screen, though, instead hav- 
ing complex filamentary and patchy structure ( Calzetti 1997 , 2001 1. 
A consequence of this is that Balmer decrement and f3 are likely to 
be sensitive to the average obscuration at different optical depths 
within a galaxy. Indeed even Ha and H/3 can be seen to probe dif- 
ferent optical depths (Serj eant et al.,|20 02). It is also the case that 
the UV spectral slope includes emission from older stellar popula- 
tions than the OB stars responsible for ionising the hydrogen (e.g., 
|Charlot & Fall|2000| |Bell|2003| >, and that these older stellar popu- 
lations are likely to be distributed differently throughout a galaxy 
than the OB stars, and hence being affected by different levels of 
dust obscuration (e.g. Calzetti 1997). In the absence of sophisti- 
cated modelling of the radiative transfer within galaxies (e.g., 



et al. 



2000 



Witt 



1992||Xilouris et al.,|1997"l|Silva et al.,|1998[|Popescu et al., 



Tuffs et al., 2004, Popescu et al.,|20TTp it is challenging to 



overcome this limitation. 

In the following sections we analyse the distribution of galaxy 
SFR, mass and redshift, on the relationship between the Balmer 
decrement, the f3 parameter and the IR and FUV luminosities. 



4 UV ATTENUATION 

We use Lir/ Lfuv as a measure of the UV attenuation under the 
assumption that the energy absorbed from the UV luminosity is re- 
emitted in the IR. We would naively expect to see higher BDs and 
steeper UV slopes with increasing Lir/ Lfuv ■ 

Figure [5] shows the variation of both Balmer decrement and 
j3 as a function of Lir/ Lfuv ■ These relationships show weak 
trends, emphasised by the best fit lines, although it is clear that there 
is a large scatter. The correlation coefficients demonstrate that these 
trends are weak at best, with r = 0.44 for Balmer decrement, and 
r = 0.48 for /3, against Lir/ Lfuv ■ It is encouraging to see that 
for low Lir/ Lfuv the values of both Balmer decrement and /3 are 
small, although at high Lir/ Lfuv the range of possible obscura- 
tion values becomes quite large (although with a clearly increasing 
mean value as Lir/ Lfuv increases). Figure [3^b) also show s the 
relationships found for local starbursts ( [Meurer et al.,|1999) Kong 



et al., 20041 and optically selected star-forming galaxies iBoissier 



et al.,|2007 ). Our sample, as also found for higher IR luminosity 



systems by |Buat et al.,] \20W\ , spans the regime between these 
models. This is obviously the origin of the discrepancies seen ( |Wi-| 
jesinghe et al., 2010, Buat et al., 2010) when applying the standard 
relations between /3 and attenuation in deriving SFRs. 

It is instructional to look at how different approaches to pa- 
rameterising the FUV attenuation, Afuv, as a function of /3 com- 
pare. We show a number of different results in Figur e[4| In addition 
to the parameterisations of Meurer et al., ( 1999) and Wijesi nghe et| 
we use our least-squares fit between Lir/ Lfuv and /3 



from Figure [3] 

Io 9 [Lir/Lfuv] =(fi+ 1.66)/0.98, (5) 
combined with Equation 2 from Buat et al., (2005), to derive the 
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Figure 3. (a) Balmer decrement and (b) /3, as a function of the IR to 
dust uncorrected FUV luminosity ratio. The thin solid lines in each panel 
show the least-squares linear best fit. Other relationships between (3 and 
Lir/ Lpuv are shown as follows. Dashed line:|Kong et al., If 2004 1 for lo- 
cal starbursts; Das h-dot line: | Me ureret a l~J^1999) also for local starbursts; 
Thick solid line: Boissier et al., (2007} for optically selected star forming 
galaxies. 



thin solid line in Figure[4] We also use the relation between Balmer 
decrement and /3 derived in §[3]above, together with the relationship 
between Balmer decrement and Afuv arising from the application 
of an obscuration curve (here we use Fischera & Dopita 2005|), to 
derive 



At 



7.49 



log ( 



0.62/3 + 6.79 
2.86 



(6) 



It is encouraging to see the consistency between the relation- 
ships of |Buat et al~"] < |2005fr and [W ijesinghe et al.,| < |2010> , and the 
offset of the relationship of Meurer et al., (1999) has been discussed 
above. It is particularly intriguing, however, to see the very different 
relationship presented when Balmer decrement is used as a proxy 
in the process of estimating attenuation. The significantly differ- 
ent slope derived here (dashed line in Figure H), a consequence of 
the logarithmic dependence on the Balmer decrement, may start to 
give some hints as to why such dramatically different results are 
obtained when using these different approaches to obscuration cor- 
rections in deriving SFRs. We return to this in §|6]below. 



5 IR AND FUV LUMINOSITIES 

Figure[5]shows a clear trend between the FUV and IR luminosities. 
The scatter within the trend varies at different FUV luminosities. 
At high FUV luminosities the distribution of IR luminosities is rel- 
atively constrained (~ 1 dex) while at low FUV luminosities the 
distribution is much broader (~ 2.5 dex). 

Figure [6] shows relationships between Lir/ Lfuv , Ha SFR, 
mass and the specific SFR (SSFR) as a function of redshift, where 
SSFR is the SFR divided by the stellar mass of the galaxy. The 
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Figure 4. Different observational correlations between j3 and the attenua- 
tion at FUV wavelength, Apuv- Thin solid line: Equation 2 from Buat et 
|aLT](2005) combined with Equation|5] Dot-dash line: |Meureretal., j!999) ; 
Thick solid line:[Wijesinghe et al., (2010) ; Dashed line: This work, Equa- 
tion[6]derived using the relationship between BD and ft. 
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Figure 5. FUV luminosity compared with IR luminosity. While there is a 
clear relationship overall between the two luminosities, at low FUV lumi- 
nosities there is a much broader range of IR luminosities compared to high 
FUV luminosities. 

range of observable values for these properties becomes more lim- 
ited with redshift, a consequence of the flux limited selection of 
the GAMA survey. This effect is strongest for SFR and mass, but 
less limiting for Lib./ Lfuv ■ The mild increase seen in average 
Lir/ Lfuv with redshift is consistent with the evolution in dust 
mass found by Dunne et al. (2010). 

In Figure [7] we compare the IR to FUV luminosity ratio, 
which indicates dust attenuation, against the IR luminosity, the 
dust uncorrected FUV luminosity, and the sum of these two quan- 
tities, in order to explore these relations as a function of SFR, mass 
and redshift. Each set of three panels in Figure [7] shows the same 
data, but colour-coded by different parameters in each row to high- 
light the impact of these properties. Presenting the combination of 



Figure 6. Variation with redshift of (a) Lm/Lpuv ratio, (b) Ho SFR, (c) 
Mass, (d) specific Ha derived SFR. The gap centred on z = 0.16 shows 
where the wavelength of the atmospheric O2 band (Fraunhofer A-line) 
overlaps with redshifted wavelength of the Ha emission line. For this rea- 
son we cannot use Ha emission line measurements that fall into this redshift 
range, leading us to omit these data from our analysis. 



Lfuv and Lir data in this way allows ease of comparison to the 
wealth of existing similar analyses (e.g., |Buat et al.,|2005| [2009) 
|20l0l >. 

It has become common now to infer SFRs by combining the 
IR and UV luminosities, as an alternative to explicit obscuration 
corre ctions (e.g.,|Buat & Xu|1996| |Flores et al.,|1999| |Buat et al.,| 
[T999"l|Hopkins & Beacom|2006> . The ranges of the values in Fig- 
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also agree with those observed in |Buat et al.,| (2009] >. |Buat] 



|et~aqp()IOl > measure the IR to FUV luminosity ratio for a sample 
of galaxies in the Lockman field surveyed by Herschel as part of 
the HerMEsj^] survey. Our sample probes lower infrared luminosi- 
ties than Buat et al., (2010), and also shows a slightly lower range 
in the ratio of IR to FUV luminosities, a consequence of the low 
redshift range of our sample. 

The basic structure seen here is that systems of low Lir have 
lower Lm/ Lfuv ratios, while for high Lir a broad range of 
Lir/ Lpuv ratios is visible, with medians increasing with increas- 
ing Lir (a trend which continues to higher Lir, as in Figure 1 of 
|Buat et al.,||2010| l. It is clear that despite the broad range of ra- 
tios at high Lir the majority of galaxies display high Lir/Lfuv 
ratios. Strikingly, the highest Lfuv systems (prior to dust correc- 
tion) are those with the smallest values of Lir/ Lfuv , a result that 
has been emphasised by |Buat et al.,| p009), and which highlights 
the limitations of UV-selection when performing a star-formation 
census. Systems with low observed Lfuv tend to have the high- 
est Lir/ Lfuv ratios, an indication of significant obscuration in 
these systems. As a function of the combined IR and FUV lumi- 
nosity, the general trend is similar to that with the Lir alone. The 
contribution of the FUV luminosity to the total is significant only 
for the lower Lir systems, and essentially negligible for the more 
luminous infrared systems. 

|Buat et al.,] ( |2009) shows that the IR to UV luminosity ratio 
does not change significantly with increasing redshift for the ma- 
jority of UV-selected galaxies. Interestingly, Lyman break galaxies 
at z < 1, which are photometrically colour selected, were found 
by Buat et al., (2009 1 to show systematically lower Lir/ Lluv at 



the same UV luminosity compared to the UV-selected systems. The 
implication here is that Lyman break galaxies maybe more biased 
against obscured systems than those identified by simple UV selec- 
tion. For IR selected (250/^m selections with the SPIRE instrument 
on board the Herschel telescope) galaxies, however, Dunne et al. 
(2010) find that the dust mass function evolves strongly with red- 
shift, so that the dustiest galaxies at higher redshifts have higher 
characteristic dust masses than those at lower redshifts. They also 
find that sub-mm selected galaxies are more dusty per unit stellar 
mass and more obscured at earlier times, however this second find- 
ing is simply based on the trend of averaged quantities over redshift 
and thus can be influenced by the IR selection which favours more 
obscured galaxies as the limiting detectable dust mass increases 
with redshift. 

Before discussing the trends with SFR, mass and redshift, we 
emphasise that this sample is not volume or luminosity limited, and 
there are consequently strong correlations between these three pa- 
rameters (with mass and SFR both being higher at higher redshift), 
simply as a consequence of the flux-limited sample being explored. 
While this selection bias should be borne in mind, none of the fol- 
lowing investigation is reliant on using a volume or luminosity lim- 
ited sample, and the conclusions are not affected by the flux-limited 
sample selection. 

The SFR dependence is shown in Figure [TJa)-(c). The SFRs 
used here are those derived from the Ha luminosity, obscuration 
corrected using the Balmer decrement and the obscuration curve 
of |Fischera & Dopita|f2005| >, as in |Wijesinghe et al.,| ( |2"0"l0| >. It is 
clear that Lir correlates with the Ha derived SFRs, and it is worth 
emphasising that systems with the highest SFRs and the highest 
Lir/ Lfuv ratios are pushed to the lowest FUV luminosities. In 
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other words, a UV-selected sample will always be prone to missing 
even very high SFR systems with sufficient obscuration. 

The redshift dependence (Figure |7jd)-(f)), and the mass de- 
pendence (Figure[7Jg)-(i)) are both very similar to the SFR depen- 
dence. This is a consequence of both the flux limit (meaning that 
high-SFR systems tend to be found at higher redshift) and the rel- 
atively tight relationship between stellar mass and SFR. We can 
account for this to some degree by looking at the dependence on 
specific star formation rates (SSFR), which are sampled more uni- 
formly with redshift (Figure |7jj)-(l)). This shows a very similar 
trend as for SFR and mass, with higher luminosity systems having 
higher SSFRs in general (a consequence of massive star forming 
galaxies being able to support proportionally higher levels of star 
formation than low mass systems). There is no strong dependence, 
though, between Lir/ Lfuv and SSFR, with a broad range of SS- 
FRs being sampled at any given value of Lir/Lfuv- 

Figure[7jd) echoes Figure 1 of jBuat et al..|pOI()i. but contin- 
ued to lower luminosity and redshift, sampling the < z < 0.35 
range. This demonstrates a clear trend that continues to lower 
Lir/ Lfuv ratios at lower Lir and lower redshift. The apparent 
strong correlation between Lir and redshift is a consequence of 
the flux-limited sample as well as the strong evolution in the Lir 
over this redshift range. Dye et al., (2010) showed that the 250/im 
luminosity density increases as (1 + z) 1 ' 1 in < z < 0.2. Sim- 
ilar evolution has also been observed from IRAS and Spitzer/ISO. 
|Dunne et al., (2010) showed that the dust content of galaxies also 
increases over this redshift range. 

It is clear, though, that lower luminosity systems, which in 
general show lower Lir/ Lfuv ratios, tend to have lower SFRs 
and masses, and are more easily visible at lower redshift. The 
tendency for higher luminosity systems (those with higher SFRs) 
to have a broader range of obscurations, and on average to be 
more heavily obscured, is now well-established (e.g., Hopkins et 



|al.|2001[|2003l|Afonso et al.,|2003||Perez-Gonzalez et al.,|2003| ) 
and this tendency is illustrated throughout Figure [7J 

The distribution of the obscuration measures /3 and Balmer 
decrement can be seen in Figure |7Jm)-(o) and Figure |7Jp)-(r) re- 
spectively. The striking feature here is the lack of any strong re- 
lationship. Reflecting the weak trends seen in Figure [3] there is a 
definite tendency to see systems with higher values of /3 at higher 
Lir /Lfuv ratios, and these systems are also those with high Lir 
and the lowest Lfuv- The Balmer decrement on the other hand 
shows even less of a systematic variation, with only a very weak 
tendency toward higher values for high Lir/ Lfuv and high Lir. 
This again is likely to be a consequence of the sensitivities of 
Balmer decrement and f3 to different optical depths within these 
galaxies. 

The overall conclusions from this exploration are: 

(i) We see an obvious relationship between obscuration as mea- 
sured by Lir/ Lfuv and SFR measured using dust corrected Ha 
luminosities. This is consistent with existing results (e.g. .[H opkins 
|et al.|2003||Afonso et al.,|2003||Perez-Gonzalez et al.,|2003| >, show- 
ing lower obscuration for low SFR systems, but a large range of 
obscuration in high SFR systems. 

(ii) We see a similar result with galaxy stellar mass, reflecting 
the mass-SFR relationship in galaxies (high SFR systems tend also 
to be high stellar mass), again with low mass systems displaying 
lower obscuration, but high mass galaxies having a broad range. 
This is consistent again with the result found by |Buat et al.,j (2009) 
as a function of K-band luminosity. 

(iii) The redshift variation is a consequence of the sample se- 
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Figure 7. The IR to FUV luminosity ratio is shown against (from left to right), the IR luminosity, dust uncorrected FUV luminosity, and the sum of the IR and 
dust uncorrected FUV luminosity. The plots are colour coded as a function of (from top to bottom) the dust corrected Hct derived SFR, redshift, galaxy stellar 
mass, dust corrected Ha derived SSFR, and Balmer decrement. (Continued below.) 



lection, with high mass, high SFR systems only being able to be (iv) The correlation between the f3 parameter and the BD was 

detected at higher redshift, due to the volume being sampled. As found to be very weak. 
|Dunne et al.,| ( [20T0| > shows, there is also likely to be a contribu- 
tion from evolution in the dust masses in galaxies over this redshift 
range. 
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Figure 7. continued. 



6 STAR FORMATION RATES 

The primary reason for a detailed investigation of the UV spectral 
slope is in assessing its utility for making obscuration corrections 
to UV luminosities, and in many cases to subsequently calculate 
SFRs. Here we look at SFRs calculated using FUV and Ha lumi- 
nosities, corrected for obscuration using various combinations of 



j3 and Balmer decrement. These SFRs were calculated using the 
prescription outlined in Wijesinghe et al., (2010) where the SFR 
conversion factor was derived using a Baldry & Glazebrook|j2003| 
IMF with a mass range of O.IMq to 120Mq and the population 
synthesis model by |Fioc & Rocca-Volmerange ( 1997) . 



10 D. B.Wijesinghe et al. 




SFR H „ (M . yr ) SFR H<1 (M.gyr- 1 ) 




Figure 8. FUV SFR plotted as a function of Ha SFRs using various obscuration corrections, (a) No obscuration corrections applied, (b) FUV SFR against Ha 
derived SFR, both obscuration corrected using the Balmer decrement, (c) FUV SFR against Ha derived SFR, both obscuration corrected using f3. (d) FUV 
SFR corrected for dust using compared against that measured using Ha, corrected using the Balmer decrement. The solid line in each panel shows equality 
between both axes. The figures show the density of points in a given region in the plot where darker the shading the denser the points in that region. 



|Wijesinghe et al.,|(2010[ l calculated SFRs from FUV luminosi- 
ties corrected for obscuration using both Balmer decrement and /3, 
showing that even when accounting for the offset in the attenuation 
from the |Meurer et al.,| ( |1999[ > prescription, there remains signifi- 
cant scatter in the relation. This is highlighted again in Figure [8] 
This Figure compares SFRs derived using dust corrections made 
using Balmer decrement for both measures, using f3 for both mea- 
sures, and using /? for the FUV and Balmer decrement for the Ha 
SFRs. The first point to emphasise is that using inconsistent obscu- 
ration measures for different SFR estimates is the primary source of 
the huge scatter seen. The scatter in Figure[8|d) is a direct reflection 
of the scatter between /3 and Balmer decrement (Figure|2](. 

It is also clear that SFRs calculated using dust obscuration 
corrections based on the Balmer decrement give the most self- 
consistent results (Figure |7jb)). When /3 is used to correct both 
FUV and Ha luminosities a correlation is seen in panel (c) in Fig- 
ure [7] This is expected as both the FUV and the Ha luminosi- 
ties share the identical /3 factor leading to the same obscuration 
factor being applied to both. Since the uncorrected luminosities 



show some agreement (Figure [8j a)), multiplying by the same fac- 
tor merely stretches out the distribution to what is observed in Fig- 
ure[5Jc). The deviation seen at high SFRs, with FUV SFRs tending 
to be slightly, but systematically, underestimated compared to Ha 
SFRs, is also visible in the uncorrected luminosities, and reflects 
the greater effect of obscuration at UV wavelengths. This contrasts 
to the Balmer decrement approach which finds agreement between 
the two SFR indicators by using dust corrections appropriate for 
each wavelength. 

Given that Lir/ ' Lfuv is (marginally) more tightly correlated 
with /3 than it is with Balmer decrement (Figure [5J, it is straight- 
forward to infer that using an infrared excess to account for ob- 
scuration at UV wavelengths will result in similar effects as seen 
above with ft. We do not attempt to calculate SFRs directly from 
Lir in the current investigation, as we are focusing on the self- 
consistency of obscuration corrections between the UV and Ha 
derived SFRs. A detailed analysis of SFR calibrations from Her- 
schel fluxes in multiple passbands is currently underway (da Cunha 
et al., in prep.). 
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Consider the implications of the scatter between j3 and Balmer 
decrement. If the origins of the scatter are largely physical, associ- 
ated with probing different optical depths, stellar populations and 
dust geometries, it is perhaps reasonable to argue that ft is a more 
appropriate metric to use in making dust obscuration corrections to 
UV luminosities, as the obscuration being probed by /? is the same 
as that affecting the luminosity to be corrected. Similarly, the argu- 
ment would be made that Balmer decrement would be the most ap- 
propriate obscuration metric to use in correcting Ha luminosities. 
Why, then does this combination give rise to a poorer comparison 
between the two estimates of SFR (FigureJSJd)) than using Balmer 
decrement for both (Figure [8jb))? We explore this further in §|7] 
below. 



7 DISCUSSION 

Given the above limitations, what can we say about the utility of 
P as a proxy for obscuration in galaxies? First, we emphasise that 
we are limiting ourselves here to values of j3 estimated from the 
two broad UV bands of GALEX. More robust constraints on the 
UV spectral slope ( jCalzetti e t al., 19941 would clearly improve the 
situation by removing that contribution to the scatter in the mea- 
surement of /3. Kong et al.. (2004) also show that the inclusion of 
SFR histories also helps to reduce the scatter when using j3 as an 
obscuration metric. The star formation rates of the current sample 
were investigated in the context of Lir/Lfuv, and also against 
/?, and showed no straightforward correlation. Regression analysis, 
using Balmer decrement as the independent variable, and multiple 
combinations of /3, SFR, redshift and stellar mass, as the dependent 
variables, shows little improvement in the weak correlation already 
found between Balmer decrement and j3. 

It is also undoubtedly the case that the UV spectral slope is 
sampling a measure of obscuration that is physically different from 
that measured by the Balmer decrement (e.g., |Calzetti|1997||Char-| 
lot & Fall 2000 1. The use of the Balmer decrement in applying ob- 
scuration corrections to both FUV and Ha derived SFRs, however, 
produces highly consistent estimates. This suggests that the effect 
of the potentially different obscurations on the UV luminosity and 
the Ha luminosity is not a major or systematic effect, and that the 
uncertainties associated with the /? approach are much larger than 
these underlying physical differences can explain. 

The conclusion we are left with, then, is that, the above limita- 
tions associated with measurement aside, the UV spectral slope is 
sensitive to many factors, of which obscuration is only one. These 
additional factors, including the age of the most recently formed 
stellar population, and contributions from older stellar populations, 
along with metallicity and IMF slope, are not insignificant, and are 
challenging to account for in a simple way. Further investigation 
of the utility of /3 as an obscuration metric will need to explore all 
these effects in a thorough and systematic fashion, ideally with data 
that samples the UV spectrum finely. 

Finally, we reiterate that the correlation of Lir/Lfuv is 
stronger with the UV spectral slope than it is with the Balmer decre- 
ment, although the correlation coefficients in both cases are low. 



8 CONCLUSION 

We have used a sample of galaxies from the H-ATLAS SDP region, 
with multiwavelength photometry and spectroscopy from GAMA 
to explore the relationships between the UV spectral slope /3, and 



the Balmer decrement. We find that there is a very poor correlation 
at best between the Balmer decrement and /3, and that the use of ft 
as an obscuration metric suffers from significant limitations. 

We see, as found by other authors, a clear but weak depen- 
dence between infrared excess, here estimated using Lir/Lfuv, 
and p. We find at most a weak trend between Lir/Lfuv and 
Balmer decrement. We also reiterate the results of Bua t et aL^I 
(20091 that UV selected samples will be strongly biased against 
heavily obscured systems, even of similarly high luminosity to 
those that enter such samples. We find consistent results with |Buat| 
et al., (2010) regarding the trend of Lir/Lfuv with Lir and red- 
shift, probing in this analysis to lower redshifts and luminosities. 
We also see trends with SFR and galaxy mass that reinforce exist- 
ing correlations and trends between these properties with luminos- 
ity and obscuration. 

In summary, we urge caution in the use of j3 as an obscura- 
tion metric when it is able to be inferred only from a few broad- 
band photometric measurements, in particular for systems at high 
redshift, given the limitations apparent in doing so even for large 
samples of well-studied galaxies at low redshift. 
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